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A system containing a pre-steady state standard (non-autocatalytic) reaction, with 
multiple paths, evolves toward a kinetic state with the minimum attainable activation 
free energy. Displacement of the path frequency distribution in this transition was 
shown to minimise the affinity linked to this change in activation free energy. In 
damping this scalar force, a standard system is driven along a path of least action, 
as previously established for a system of competing autocatalytic reactions. A kinetic 
source of time asymmetry arises within the system, as the activation affinity moves 
the system toward the most probable distribution of reaction paths. As the functions 
of state are not changed by path displacement, a change of kinetic state cannot 
produce chemical work. This generalises the notion of force to a scalar quantity 
responsible for a displacement that does not yield work or heat. Spectrophotometric 
observations reported on the transition to steady state kinetics during dinitrophenyl 
phosphate phosphorolysis confirmed that time variations in the activation affinity are 
non-positive. 
W 
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1 . Introduction 

A deviation from tine mean activation free energy for nucleotide condensation in a 
lieterogeneous population of competitively replicating polymers results in a path- 
centred force (Davis, 1994, 1996) that shifts the polymer frequency distribution and 
increases the rate of synthesis. Elimination of polymer species with suboptimal 
template competence accounts for this increase (Kramer et al., 1974; Eigen et al., 
1989). In effect, nucleotide condensation tends to take place by progressively more 
rapid reaction paths, or, equivalently, by more probable transition states. The system 
thus evolves with time in the direction of the most probable kinetic state. 

This leads to the proposition that two independent sets of forces can orient a 
chemical system in time. A scalar force that increases internal entropy accounts for 
the thermodynamic orientation of a process in time. The time asymmetry in 
competitive replication, and putatively in other pre-steady state chemical systems 
(Laidler, 1987), however, arises from the spontaneous displacement of path 
frequencies between a reactant and product state, which remain unaltered. In 
addition to thermodynamic forces responsible for changes of state, therefore, a 
chemical system also apparently contains path-centred forces that move it toward 
the most probable set of paths between states. 

It appeared of interest to establish that the proposed theory can be extended 
to a change of kinetic state by a standard reaction. This would reveal whether 
changes in a standard pre-steady state system (i) damp the activation affinity at a 
maximal rate, and (ii) yield a path of least action, as in a set of autocatalytic 
reactions. In contrast to evolution among competing autocatalytic reactions, where, 
in principle, only the fastest reaction path is ultimately retained (Davis, 1996), all 
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available paths are retained during the transition to steady state kinetics by a 
standard reaction. Moreover, steady state kinetics are achieved rapidly in the latter. 
A transition interval of less than 200 ms, for instance, occurred in an enzyme 
catalysed reaction (Trentham & Gutfreund, 1968) considered here. This compares 
with a transition interval of nearly 3 hr in the RNA evolution experiment (Kramer et 
al-, 1974) investigated previously. For the same principles to govern a change of 
kinetic state in two systems with such significant qualitative and quantitative 
differences, they would have to apply in general. 

2. Transition to Steady State and the Activation Affinity 

A standard reaction with more than one transition state, T(i), forms a set of parallel, 
consecutive reactions, 

k k 

A ^ {TJ ^ B ; i = 1,2,...,n; k^ » k„ k, >> k^ (1) 
k-1 k.2 

Each reaction path is seen to produce and consume an unstable intermediate, with 
rate coefficients k|. The frequency distribution of the set {T(i)} determines the 'kinetic 
state' of system. For definiteness, a closed, isothermal system at constant volume is 
considered. 

A change of kinetic state within the system may be linked to the presence of a 
scalar force - the activation affinity. It reflects the difference in activation free energy, 
F*, of a reaction path at time, t, from its stationary value, some time later (Davis, 
1996): 
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A*(i) = -(F*33(i) - F* (i)) = RT(ln K^i) - In K*,(i)) 

= -RT In f*(i) (2) 

R is the gas constant and T temperature. K* is an activation state 'equilibrium 
constant'; K*(i) = T(i)/A. A time-dependent expression relating the i*^ transition state 
concentration, Tt(i), to its steady state level, Tss(i), is obtained on solving the rate 
equations for (1); 



f*(i) = K* (i)/K33(i) = T/T33 

= k2j_(1 -e-(''''-'^''") (3) 



Aq is the initial reactant concentration. Steady state requires i » leaving the 
factor (1 - e"*'^'^ ■''^'''') as a measure, f*(i), of advancement by the i'^ reaction path to 
steady state. When the reaction (1 ) is catalysed 

A + C ^ {T(i)} ^ B + C 

K*(i) is T,(i)/(AC) and K33*(i) equals T33(i)/(AC), which reduces to k^lk^. The second 
order advancement term, f*(i)'^', is 
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(k2,-k,,(A+C))' 



(1 - J<2,i^i,i(A + C}_ exp-[k2,it - H dt k,,i(A + C)] 
\ k2,i - k, i(Ao + Co) 



where kgj » k^ i(A + C). 



The steady state constraint, requiring that the first step in (1) be rate limiting, is a 
feature of the reaction. It implies that the formation of each T(i) involves reactant 
activation. In the second step of (1), T(i) decomposes rapidly to product, at a rate 
determined by 'diffusion' along the high energy segment, 6, of the reaction path. The 
decomposition frequency of the i'^ transition state is 



\2nmi*/ 6i 

kgT being the thermal energy per molecule, and m*| is the mass of the i'^ activated 
complex. Although kg, and k^ | appear equal in (1), reaction rate theory (Eyring and 
Eyring, 1963) specifies that, once formed, the activated complex proceeds to 
product. A transmission coefficient less than one being discounted here, and 
quantum tunnelling through the activation energy barrier is disregarded also. When 
the reaction (1) proceeds spontaneously in its forward direction, each formation rate 
coefficient, k^ exceeds the corresponding reverse coefficient, k^,. The former being 
defined here as. 



h is Planck's constant. Q*,, are the partition functions of the i'" transition state 
(omitting the translation term; (2nm|*kBT)'''Vh) and reactant; QaQc replaces for the 
catalysed reaction path. i denotes the zero point energy per mole. In addition to 




ki,i = (2iimi!k,BWe-' 
h Q, 
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the free energy of activation, implicit in tlie relation for |, K2,i reflects the free 
energy difference between A and B. 

The rate of a reaction, in which all paths are retained during a change of kinetic 
state, is increased by addition of a path. Path number ordinarily increases with the 
introduction of a catalyst. Total transition state concentration increases on passing 
from an uncatalysed, T^j, to catalysed, T^ + T^,, kinetic state. The change in activation 
free energy is 

AF* = F*, - F*, = -RT In [1 + KVK*J < (4) 

Kj,* = JJ{AC) and K^j* = T^A, where T, A and C are transition state, reactant and 
catalyst concentration - activity for (4) non-ideal. The inequality applies for all values 
of the activated 'equilibrium constant' K*; K*^, I K*^. A decrease in overall activation 
free energy requires only that the number of paths increase; transitions with a large 
number of paths are conceivable. Conversely, a spontaneous decrease in path 
number is excluded by (4). This would require that a standard system 
spontaneously shift to a less probable set of paths. The probability of this transition 
decreases exponentially with the increase in activation free energy. A spontaneous 
transition among competitive autocatalytic processes, by contrast, will exhibit a 
significant reduction in path number. As slow paths are eliminated, the most 
probable kinetic state then consists solely of the fastest reaction path! 

A spontaneous change of kinetic state is characterised by a positive activation 
affinity. From (2) and (3), 



e-A*(i)/RT ^ T,(i)/T33(i) < 1 , A*(i) > 



(5) 
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A positive kinetic force implies a transition state lias a probability in the canonical 
distribution of states that is higher than its pre-steady state probability. With 
relaxation to the canonical distribution during reaction advancement, transition state 
probability increases to its steady state level and the kinetic force vanishes. A 
positive force of path displacement does not represent a potential for chemical work, 
since system free energy is not altered by changes in path probability between A 
and B states (1). A change of kinetic state thus involves a novel kind of force: a 
scalar quantity linked to a displacement, which produces no work or heat. 

3. Force Damping by Pre-Steady State Changes 

It can be demonstrated that changes in path probability in a pre-steady state system 
optimally damp the activation affinity. 

Differentiation of (2) with respect to time reveals that the activation affinity is 
damped as 

A*(i)' = -RTcl)i < 0, > (6) 

cl)i is (k2,i - ki i)e ''^'■'-'^■'"/(l - e f^^^ ^Li)*) in a first order reaction and 

(k^, - K-iA + O) exp-[k2,t - dt k,/A + C)]/(1 - exp-[k,,t - dt k,/A + C)] 

at second order. As noted earlier (section 2), k2i>> k, , in the former and kg, » k^ |(A + 
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C) in the latter is a precondition for acliieving steady state kinetics. In view of (2), the 
damping rate (6) corresponds to d(ln f*(i))/dt. Since f*(i) reduces to T,(i)/Tg3(i), as T(i) 
approaches a steady state concentration, Tss(i), both the activation affinity (2) and its 
damping rate (6) diminish and ultimately vanish. 

The probability that a reaction will occur by its i'" path at time t is 

p,(i) = K*(i)/EiK*(i) (7) 

Its rate of change is 

p,(i)' = p,(i)cDi(1 - Si) (8) 

In this expression for Pt(i)', s, is a control variable, whose value determines the 
direction of change in path probability with time. S| equals where $ is 
E| ct)iP,(i). It is demonstrated in Appendix A, by the method of undetermined 
multipliers, that this rate equation maximises $ for a path frequency distribution of 
given variance. It follows from (6) that the mean activation affinity is optimally 
damped. 



4. Least Action Path to Steady State 

The action specified by mean activation affinity values (2) over time in a pre-steady 
state standard reaction is 
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S = rodtA* 

= -RTj; dt p,(i) In [f*(i)] (9) 

Path frequency changes (8) accompanying the transition to steady state are 
demonstrated in Appendix B to give a stationary value for the action. It was shown 
above (section 3) that the activation affinity is minimised at each instant during this 
transition. It is concluded from this finding that a stationary value for the action 
(Appendix B) corresponds to the system following a path of least cost in terms of 
action, while evolving toward steady state. 



5. Analysis of a Pre-Steady State Reaction 

Spectrophotometric measurements of 2,4-dinitrophenol (DNP) formation on 
phosphorolysis of 2,4-dinitrophenyl phosphate (DNPP), in the presence of alkaline 
phosphatase (EC 3.1.3.1), using the stop-flow technique (Trentham and Gutfreund, 
1968) illustrate the effect of a change of kinetic state on the kinetic force in a 
standard reaction. 

Figure 1 shows the rate of DNP formation at 25°C and pH 5.9 decreased over 
time, with an effective first order rate coefficient of nearly 30 s"". A decreased rate of 
DNP formation accompanied reaction advancement, because dephosphorylation of 
the enzyme-phosphate complex, E-P, 
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k k 
DNPP + E ^ DNPP.E - E-P + DNP 

K 

E-P + H2O - E + Pi (10) 

proved to be rate limiting under acidic conditions: kg « k^ or kg. Consequently, E-P 
levels increased at the expense of active enzyme molecules, E. Uncatalysed 
hydrolysis of the phosphate bond in DNPP 
K 

DNPP + HgO^ DNP + Pi (11) 

to produce DNP and an inorganic phosphate ion. Pi, did not exceed 1 per cent of the 
catalysed rate. 

FIGURE 1 

As the transition state in the catalysed reaction, E-P.HjO, increased toward its 
steady state concentration, the activation affinity, A*^, decreased. From a value of 
nearly 14 kJ mol"" at 0.1 ms, the force reduced to only 0.03 kj mol"' by 150 ms. Its 
rate of decrease (6) 

AV = -RT30e-^°V(1 - e"^°') (12) 

also diminished with time. The rate was -2440 kJ mol"* s"' at 1 ms. At 10 and 100 ms 
after the start of the reaction, it reduced to -213 and -4 kJ mol"' s"\ respectively. With 
the uncatalysed reaction (1 1 ) remaining at its steady state, the activation affinity for 
this path, A*^j, was zero throughout the experiment (Fig. 1). 

The action based on mean activation affinity (9) values during DNPP 



phosphorolysis was evaluated. The probability of catalysed and uncatalysed paths 
for the pre-steady state interval relied on a steady state probability for the latter of 
0.01 . An interval of 153 ms was required to reach a catalysed path probability close 
to 0.99. The action over this interval was 87 kJs rr\o\'\ 

Results presented earlier (section 2,3) establish that activation affinity damping 
rates and the action obtained for the DNPP reaction, in its transition to steady state 
kinetics, are extremal values. 



6. Discussion 

Changes in the path probability distribution of a standard reaction, with multiple 
paths, could have been investigated at equilibrium. In this state, no thermodynamic 
force acts on the system and no change occurs in reactant and product 
concentration. The system is then effectively in a state of maximum symmetry 
(Callen, 1974, 1985). In particular, it exhibits no time asymmetry of thermodynamic 
origin. Addition of a catalyst to a chemical system at equilibrium, notwithstanding 
thermodynamic invariance, elevates the rate of interconversion between reactant 
and product molecules, to produce an effect similar to that observed on mixing 
alkaline phosphatase and DNPP (section 5). In the presence of a catalyst the 
system moves to a more probable kinetic state, as the path probability distribution of 
a reaction shifts to lower the free energy of activation (section 2). 

During this transition, a chemical system at equilibrium changes in a preferred 
direction with time. Since this time asymmetry arises in the absence of 
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thermodynamic forces, a second source of time asymmetry evidently occurs witfiin 
the system. A second scalar force was identified: the activation affinity. This finding 
adds to the number of known sources (Davies, 1974; Penrose, 1979; Schulman, 
1996) of time asymmetry. 

Both thermodynamic and kinetic time asymmetries derive from statistical 
displacement of a system. They reflect a probability difference between an initial and 
final state, and forces of both types are damped as the system moves toward its 
most probable state or set of reaction paths. At the most probable state consistent 
with kinetic and thermodynamic constraints, these forces vanish and no further 
change occurs. The amount of heat contained by a system, at a specified 
temperature, is then maximised, in accord with the second law of thermodynamics. 
In addition, kinetic theory requires that the rate of molecular interconversion is also 
maximised, for a specified set of reaction paths. 

The activation affinity displaces a chemical system, without yielding either work 
or heat. Shifting the path probability distribution between reactant and product state 
causes no change in system free energy and, therefore, no change in work potential 
(section 2). The concept of a kinetic force thus generalises the notion of force. It 
requires inclusion of a scalar quantity responsible for a kind of displacement, which 
is not coupled to work output. Similar to classical and quantum forces in a 
conservative dynamical system (Feynman & Hibbs, 1965; Soper, 1976) and the 
thermodynamic force (Shiner & Sieniutycz, 1994), the kinetic force was 
demonstrated to move a system along the path of least action (section 3). 

Danwinian evolution based on natural selection among competitively replicating 
RNA molecules (Kramer et al., 1974) was shown to damp the activation affinity for 
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nucleotide condensation (Davis, 1996). Tliis system evolved toward steady state 
kinetics over an interval of 3 hr. In this time, the effective activation free energy for 
nucleotide condensation decreased from 86.90 to 86.55 kj mo\'\ Consistent with the 
time dependence of the mean polymer formation rate coefficient (Davis, 1978), the 
mean activation affinity, varied with the length of the replication interval, L, 

A* = 0.33 + 3.27x10"^ L - 1.35x10"^ + 7.15x10"^^ 

decreasing from 0.33 to 0.0076 kJ mol'^ over 190 min. (L = 1 1400 s). The chemical 
affinity driving nucleotide condensation during polymerisation was held virtually 
constant at 32.2 kJ mor\ throughout the transition to steady state kinetics. 
Danwinian evolution under defined physical conditions thus damped the kinetic force 
within a chemical system, moving the system toward the most probable (fittest) 
reaction path. The principle of natural selection is seen to complement Boltzmann's 
evolution principle, requiring that state-centred forces are damped with time by a 
statistical mechanism. 
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Appendix 

A. Maximum Rate of Increase in Transition State Concentration 

Changes in reaction patli frequency witli time in a pre-steady state standard system 
shall be demonstrated to maximise the rate at which the concentration of each 
transition state approaches its steady state level. 

Recalling (7) and (8), path probability and its rate of change are 

p,(i) = K* (i)/Ei K* (i) 

p,(i)' = p,(i)cDi(1 - Si), Si<1 - p,(iy>0 (A1) 

Si > 1 - p,(i)' < 

As indicated, s, controls the direction of change in path probability. K* is the 
transition state to reactant concentration ratio at time t, and (^^ is d In K* (i)/K*ss(i)/dt. 
The latter relates changes in the current level of the i'^ transition state to its steady 
state value. Should path frequency changes (A1) maximise $, pre-steady state 
changes in a standard chemical system would advance transition state levels toward 
steady state at their maximum rate. This implies maximum damping of the activation 
affinity (6). 

Displacement of the path frequency distribution is subject to two constraints: (i) 
the probabilities are normalised, so E| 6p(i) = 0, and (ii) the variance in cf) is 
momentarily fixed, E| p(i)(ct)| - = B, where B is a constant. Using (A1), constraint 
(ii) can be recast as 



Ei 6p(i) 5ln p(i) = B 
6t 6t 
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We may then define the function 



4J = Ei cl)i6p(i) + Ai(Ei 6p(i) 5ln p(i) - B6t) + K^L-, 6p(i) 



(A2) 



5t 6p(i) 



It is stationary with respect to variations in probability, 6p(i), when 



aijj = + 5p(l) + A2 = 
55p(i) p(i)6t 



(A3) 



The Lagrange multipliers are found to be = -V2 and = Substitution into (A3) 
reveals 



The probability variation is a maximum, as differentiation of (A2) again with respect 
to 6p(i) reveals the second derivative is negative: 

aV55p(i)' = - l/p(i)6t < 0, < p(i) < 1 and 6t > 0. 
B. Transition Path of l\1inimum Action 

An action (cost function) for the pre-steady state transition by a standard reaction 
with multiple paths (1) is given by 



(6p(i)/6t),,, = p(i)(cDi - $) 



(A4) 



S = - RT J^odt p,(i) In [f*(i)] 



(81) 
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the integral of the mean activation affinity (2) over the transition interval. During this 
interval, the path probability distribution (Pt(1),p,(2),...) approaches its steady state 
(Boltzmann) distribution. Using Pontryagin's minimum principle (PMP) from control 
theory (Pontryagin et al-, 1962), changes in path probability specified by the state 
equation (Al) will be shown to yield a stationary value for the action (B1). 
The cost variable adjoined to the state equation (A1 ) is 

p,(0)' = - RT p,(i) In [f*(i)] (B2) 

The resulting Hamiltonian from (Al) and (B1) is 

H = - RT K,E, p.(i) In [f*(i)] + p,(i)cDi(1 - s-) (B3) 

The Hamiltonian relations then yield 

Ao' = - aH/ap,(0) = - Ao = 1 (PMP) 

A- = - aH/ap,(i) = RT In [f*(i)] - K,<^,0 - S;) (B4) 



Aq is set equal to 1 on applying PMP, after Hamilton's relations show it is a constant. 
Solving the ODE in (B4) for A| time dependent yields 

h,{t) = exp [-Ji <\>,{^ - s-m Jo RT In [f*(i)] exp [Ji <\>,0 - sm dt 

Substitution for A| in (B3), for all i, gives a time dependent Hamiltonian. Since it is not 
constant for the pre-steady state transition path of the system, it is discounted as a 
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non-optimal path (Pontryagin et al., 1962). 



For A| independent of time {K- = 0), (B4) gives (Davis, 1996) 



Aj = RT In [f*(i)] 
cDi(1 - Si) 



(B5) 



Both RT In [f*(i)/(|)i(1 - S;) and its time variation (d(RT In [f*(i)])/dt/d((t)i(1 - Si)/dt) are 
equal to some constant, c. Recalling the relation between (^^{^ - s) and the rate of 
change in path probability, p(i)', in (Al) and the definition of the kinetic force (2), it is 
apparent from (85) that 



Here, the probability of a reaction path at some time, t, during the transition to 
steady state depends on its initial probability, Po(i), and the applied kinetic force, 
A*|(t), over the transition interval (0,t). The constant c has been equated with -RT, 
following inspection of (B5) and on noting that f*(i) measures advancement of the i 
reaction path toward steady state or, equivalently, changes in probability of the i'^ 
reaction path. Substitution of the time invariant costate variable, A|, in the 
Hamiltonian (B3) reveals 



c d In p(i)/dt = - A* 



(B6) 



With integration 



p,(i) = Po(i) exp [J^ A*(t)/RT dt] 



(B7) 
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H/RT = - p,(i) In [f*(i)] + p(1 ) In [f*(1 )] + p(2) In [f*(2)] + . . . (B6) 
H = 



The control variable, Si, specified by (B5), therefore, yields a stationary Hamiltonian. 
Control theory indicates that a zero-valued Hamiltonian, over an interval with a freely 
chosen final time, t > 0, corresponds to a stationary value of the cost function. As 
(B1) is the cost function for the pre-steady state transition, the path traced yields a 
stationary action. It was previously established (Appendix A) that the activation 
affinity, on which this action (B1) is based, is minimised at each instant throughout 
the transition interval. Consequently, a pre-steady state system can be concluded to 
follow a path of least action. 
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Figure 1 . Damping of the activation affinity, A*, by the shift to steady state kinetics 
during phosphorolysis of DNPP. Reaction advancement was determined by 
spectrophotometric measurements (Trentham and Gutfreund, 1968) on DNP 
liberation over an interval of 175 ms, after initiation of DNPP hydrolysis by alkaline 
phosphatase at 25°C and pH 5.9 with the stop-flow technique. Continuous line, 
enzyme catalysed. Broken line, uncatalysed. Points, transition state concentration 
as a fraction, f*, of its steady state concentration. These values were calculated for 
indicated times, t, on the spectrophotometric record of DNP release. The f* values 
obtained correlate with the time course of a first order process (continuous line) 
having a rate coefficient of SOs"' (r^ = 0.976). 



